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Ph! Abstract 

Oh: 



We discuss constraints imposed on the zero-recoil curvature and higher deriva- 
tives of the Isgur-Wise function by a general quark model. These constraints are 
expressed as bounds for a given slope parameter, and compared with those based 
upon analyticity properties of QCD spectral functions. Our results also indicate 
that in the analysis of the experimental data for semileptonic B —>■ D^*> decays it 
may be important to include at least the third term in the form factor expansion 
about the zero recoil point. 



1 Introduction 

It has been widely recognized for some time that heavy quark symmetry (HQS) |l|, 0| 
enormously simplifies the analysis of the semileptonic B -^ D^*^ decays. The six form 
factors needed for the description of these decays are in the heavy quark limit reduced to a 
single unknown form factor, the Isgur-Wise (IW) function C,{uj), where u = f^ -f/jc.) is the 
product of the four- velocities of the two mesons. Furthermore, HQS also provides us with 
a prediction for the normalization of the universal form factor at the zero-recoil point, 
i.e. ^(1) = 1. As a consequence of that, normalizations of the physical form factors G{u!) 
(for B —>■ D decay) and J-'{uj) (for B —>■ D* decay) ,[] are determined up to radiative and 
power corrections. Therefore, by extrapolating the experimental data for the differential 
decay rate to cj = 1 one can obtain an accurate measurement of the Cabibbo-Kobayashi- 
Maskawa parameter \Vcb\- The decay B -^ D*lui is ideally suited for this purpose [Q. It 
is experimentally clean mode, and the decay rate at zero recoil is is protected by Luke's 
theorem against first order l/mg corrections 0]. 

This analysis has already been performed by several experimental groups 0-0. In 
IP and the fit to the data assumed a linear form for J-'iuj), i.e. J^{uj) = J-'{1)[1 — 
a{uj — 1)]J3 while in and f^ fits with quadratic form of the hadronic form factor were 
also attempted, but with the conclusion that with the existing data samples it was not 
possible to distinguish between the linear parametrization and those with more degrees 
of freedom. It should be obvious that retaining only the first term in the expansion will 
inevitably lead to an underestimate of the slope parameter. This point was made some 

time ago by Burdman [Q , who included the curvature (quadratic) term in a two parameter 
^In the absence of symnietry breaking corrections these form factors would coincide with the IW 

function ^(w). 

^In PI several other paranietrizations for J-'{uj) were also used, but all with only one degree of freedom. 



analysis. Because of the statistical uncertainty introduced in a two parameter fit, it is 
clearly important to obtain some theoretical insight about the expansion parameters, in 
order to guide the extrapolation to cj = 1. 

This issue has already been addressed in several papers [rU|-[r^|, by employing analyt- 



icity properties of QCD spectral functions and unitarity. The resulting bounds proved to 
be weak due to the presence of the T poles below the BB threshold (or possible B^ states 
below the BD^*^ threshold in [0). In the most recent work [Q Caprini and Neubert 
(CN) improved bounds for the zero-recoil slope and curvature (i.e. the second term in the 
expansion) of J-'{uj) and Q{uj) by identifying a specific B ^ D form factor which does not 
receive contributions from the ground state B^ poles. These authors have derived con- 
straints between the slope and curvature of that form factor using analyticity properties 
of QCD spectral functions and unitarity, and then used heavy quark symmetry to relate 
these results to corresponding constraints for J-'{uj) and Q{uj). 

In this letter we discuss an alternative approach for obtaining allowed regions for the 
curvature and higher derivatives of the IW function ^(cj) for a given value of the slope. 
Our results are obtained in the heavy quark limit and in the valence quark approximation, 
with some physical input about the shape of the heavy-light meson wave function. An 
advantage of the method is that, given the above assumptions, constraints can be obtained 
not only for the curvature, but also for any higher terms in the expansion. Even though 
we do not take symmetry breaking corrections into account, we believe that our results 
may also shed some light for guiding the experimental extrapolation to uj = 1, especially 
for estimating the higher order terms in the form factor expansion. 

The rest of the paper is organized as follows: in Section ^ we present the general 
valence quark model expressions for the IW function, and also for the particular terms in 
its expansion about uj = 1. In Section |^ we show how to extract bounds on the higher 
expansion parameters, if the slope of the IW function is given. Results are discussed in 



Section ^, and conclusions are presented in Section ^. In particular, we conclude that in 
the case oi B ^ D^*^ semileptonic decays an expansion of the IW function about the zero 
recoil point is converging slowly for u close to the maximum velocity transfer cUmax-Q 

2 IW form factor in a general quark model 

In the valence quark approximation the expression for the IW function describing the 
B -^ D^*^ transitions, is given in terms of the S-wave radial wave function R(r) and 
energy E of the light degrees of freedom in the ground state heavy-light meson [p!5l-||19||. 



eM = -^(Mkr)) , (1) 



where jo is the spherical Bessel function and 



The expectation value {F{r)) is given by 

(F(r)) = / r'^dr\R{r)\'^F{r) . (3) 
Jo 

We define expansion of ,^(ti;) around tu = 1 as 

eH = l-a{uj-l) + b{uj - 1)2 - c{uj - If + d{u-iy + ... . (4) 

Using (|^) it is straightforward to find expressions for the slope a and higher order terms 
in (Q) [l7, IS]. We list here the first four terms: 

a = l + l^'ir'), (5) 

b = l + lE'{r') + ^E'{r'), (6) 

c = l + lE2(r2) + lE^(r^) + ^E6(r6), (7) 

8 4 ^ ^ 20 ^ ^ 630 ^ ^ ^ ' 

d = L + -EHr^) + —E\r^) + —EHr^) + -^—E\r^) . (8) 

16 6 ^ ^ 20 ^ ^ 315 ^ ^ 22680 ^ ' ^ ' 



^For B — > D*lvi ujmax — 1-5, and for B -^ Dlvi decays uJmax — 1-6. 



Note that all quantities are positive definite. From (^ one finds E^{r'^) = 3(a— |), which 
can be used to reexpress (|), (^ and (j^) in terms of a as 

c = - + -(a--) + —f3(a--y + —^(a--)\ (10) 

8 4^ 2^ 20^^ 2^ 70 '^ 2^ ' ^ ' 

d = — + -{a--) + —p{a--f + —^i{a~-f + —5{a--f. (11) 

16 2^ 2^ 20^^ 2^ 35'^ 2' 280 ^ 2^ ^ ^ 



Here, we defined dimensionless quantities 






P = hy.^ (12) 

7 = 7::^, (13) 



From (^ one can see that a > 1/2, and therefore it is immediately evident that all of the 
above parameters must be positive.^ However, we can bound them more stringently in 
order to yield more useful restrictions on the allowed ranges of b, c, etc. 

3 Bounds 

Since all terms which accompany /3, 7 and 6 in (P)-([TT|) are positive definite, it should 
be obvious that by restricting the allowed range for those parameters we also restrict the 
range of allowed values of b, c and d, for a given slope parameter a. In other words, we 



^The bound a > 1/2 is a consequence of the prefactor 2/(w + 1) in d^), derivation of which is discussed 
in depth in |19|] , and which is closely related to the valence quark approximation. The relationship between 
the HQET sum rules and quark models was investigated in |g^, where it was shown that Bjorken ||2^, ^ 
and Voloshin |23] sum rules can be used to construct a model which is self-consistent in the heavy-quark 
limit. 



want to find PminiPmax) so that 

Pmin J^ P —1 Pmax ) V -'-'-'/ 

and similarly for other parameters. 

Without making any further assumptions about the particular form of the heavy-light 
wave function, lower bounds for (3, 7 and 6 can be obtained by considering the Schwartz- 
type inequality 

(^2m^^2n _ ^rYf) > ■ (16) 

For m = 0,1 and n = 1,2 this yields 

Pmin -L ; V ^ ' / 

7mm = 1 , (18) 

Srmn = 1 • (19) 

In order to estimate the upper bounds we need some physical input. Let us for the 
moment assume that the ground state wave function of the light degrees of freedom in a 
heavy-light meson is given in the form 

R{r) oc exp (-r'^) , (20) 

where k > 0. Note that any scale or normalization dependence in the wave function is 
unimportant, since it would cancel out in the ratios (|T2|)-(|l^. For example, k = 2 would 
correspond to the harmonic oscillator wave function, which is (with appropriate scale 
parameter) often used as an approximation for the meson wave function |2^, ^. Case 



k = 1 corresponds to the pure exponential, which seems to be favored by lattice QCD 
]. Using ( pOD one can find the expression 



which with n = 2, 3 and 4 yields j3, 7 and 5 for any k. In Table |l| we list the actual 
numbers for several cases of interest. Clearly, as the value of k gets smaller, parameters 
/3, 7, and 5 get larger. Therefore, the smallest acceptable value of k will lead to the 
largest acceptable values for our parameters. Since lattice simulations ||2^ support a pure 
exponential falloff of the meson wave function (/c = 1), one might argue that choosing, 
for example, kmin = 1/2 would leave more than enough room for possible uncertainties in 
the specific choice (|^) of the long distance behavior of the wave function. In that case 
we would have (see Table |l]) 

/Smax ^ 5.67, (22) 

^ma. ~ 107.19, (23) 

5^a. ^ 5091.38. (24) 

We wish to emphasize here that any value fcmm < 1 would be acceptable as far as this 
part of the analysis is concerned, because it is essentially guided only by the information 
obtained from the lattice [^. By choosing kmin closer to 1, one would obtain quite narrow 
range for all parameters under consideration, as we shall see in the following section. 

4 Results 

Let us first discuss the second term in the expansion (^). In Figure Q with full lines 
we show the acceptable range for the curvature 6 as a function of the slope a, with a 
particular choice of kmin = 1/2. The lower bound (denoted by L.B.) follows from (jl^), 
and the physically motivated result with fc = 1 is shown as well. With dashed lines we 



further show the result of the analysis performed in [14|, which was also obtained in the 



heavy quark limit, but (unlike ours) includes short- distance corrections. To be completely 



clear, we give here (in our notation) the form factor definition used in \T^ 



f (cu) = 1(1) [1 - d{iu - 1) + biw - 1)2 + . . .] . (25) 

To avoid confusion, we have used tilde with the CN expansion parameters and form factor. 
When the short- distance corrections (^(1) — 1.02) are neglected, (^) coincides with (§).0 
Since the CN ellipse shown in Figure |l] is rather narrow, these authors also give the 
approximate relation between the slope and the curvature as 

b - 0.72a - 0.09 . (26) 

On the other hand, our result with k = 1, yields 

6 ^-0.06 + 0.25a + 0.75a2 , (27) 

which is within the CN bounds for values of a smaller than about 0.7, and grows faster 

than ( p6D with increasing a. From the Figure |1| one can see that given a value for the 

IW function slope a, the valence quark model yields a range for the curvature b which is 

comparable in size to the CN bounds, but with somewhat higher values for b when a is 

greater than about 0.7. We remind the reader that (|) requires a > 0.5. 

The CN approach is expected to break down for higher than the second terms in the 

expansion of the IW form factor. The reason is the possible presence of the sub-threshold 

singularities which are due to scalar Be resonances, or due to states of the form {B^*^ + h), 

where h is a light hadron. We can, however, estimate the acceptable range for those 

terms in the same way as we did for the second term b. Results for the third and the 

fourth term (c and d, respectively), are shown in Figures ^ and |^. Naturally, for higher 

expansion parameters the uncertainty is increasing. Nevertheless, if kmin were close to 1, 
^We note here that perturbative corrections are expected to be at most 10-15%. For J-{uj) in [Q it 
was found that a ::i a — 0.06 and b :^ c — 0.06 — 0.06a, while for G{iu) corresponding results were found to 
hea~a + 0.02 and 6 ~ 5 + 0.01 + 0.025. 

8 



the range of acceptable values for c and d would be quite narrow. We note results for c 
and d obtained for the physically motivated case of A; = 1, where we find 

c ~ -0.03 + 0.380^ + 0.50a^ , (28) 

d - -0.01 - 0.03a + 0.09a2 + 0.38a3 + 0.31a^ . (29) 

One should also observe that all expansion parameters are roughly of the same order of 
magnitude, so that the only suppression for n-th order term in the expansion (^ is due 
to a factor of {u — l)". This fact is best illustrated by taking a ~ 1 in (^7]), (|2^) and (^9]) , 
which leads to 6 ~ 0.94, c ~ 0.85 and d ~ 0.74. Taking these values near the maximum 
velocity transfer in B ^ D*Wi decays {umax — 1-5), we find from (0) 

^(1.5) ~ 1 - 0.50 + 0.235 - 0.106 + 0.046 - . . . . (30) 

Although the curvature term is large (~ 0.235), the subsequent terms are not negligible, 
and it is obvious that series converges slowly for u ~ ujmax- 

In order to show the effects of increasing of the number of terms in the form factor 
expansion about u; = 1, in Figure ^ we show what happens as we include one, two, 
three, and four terms in (^, for A; = 1 and a = 1. Clearly, in this case keeping only 
the first two terms (a and b non-zero) in (^ is an excellent approximation to ^{uj) for 
uj < 1.2. However, as u increases, the higher order terms make a difference as far as 
the shape of the form factor is concerned. For slope values smaller than one the form 
factor convergence is better. This is illustrated in Figure |^, where we have used a slope of 
a = 0.75. Nevertheless, in the analysis of the experimental data one has to keep in mind 
that for larger values of the slope it may still be important to include at least the third 
term in the form factor expansion about u = 1. 



5 Conclusion 

Within the framework of the general quark model, we have addressed the issue of the 
bounds on curvature and higher derivatives in the expansion of the Isgur-Wise function 
about the zero recoil point. These terms are important in experimental extrapolation of 
the form factor towards u = 1. Except for slopes in the range of 0.5 to 0.7, our results 



indicate slightly larger curvature than the one obtained by Caprini and Neubert [|I4[. We 
also find that including a third term in the form factor expansion about uj = 1 may be 
important in the analysis of the experimental data. 
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TABLES 



Table 1: Parameters {3, 7 and 5 for /c = 2, A; = 1 and k = 1/2. 



k 


P 


7 


5 


2 


f (^ 1.67) 


f(^3.89) 


f (^ 11.67) 


1 


|(= 2.50) 


f (^ 11.67) 


ip(= 87.50) 


1/2 


|i(- 5.67) 


^(^ 107.19) 


^Yfg^f (^5091.38) 
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FIGURES 

Figure 1: The valence quark model prediction for bounds on the IW function curvature b 
in terms of its slope a (full lines). The hne denoted by kmin = 1/2 represents the upper, 
while the line denoted by L.B. represents the lower bounds. The dashed line shows the 
result of the CN analysis 



Figure 2: The valence quark model predictions for bounds on the third term c in the 
expansion of the IW function (P, in terms of the slope parameter a. The line denoted 
by kmin = 1/2 represents the upper, while the line denoted by L.B. represents the lower 
bounds. 

Figure 3: The valence quark model predictions for bounds on the fourth term d in the 
expansion of the IW function (p, in terms of the slope parameter a. The line denoted 
by kmin = 1/2 represents the upper, while the line denoted by L.B. represents the lower 
bounds. 

Figure 4: Illustration of the effects of increasing the number of terms in the form factor 
expansion about u = 1. We used a = 1 and the physically motivated case of /c = 1. 

Figure 5: Same as in Figure H except with a slope parameter a = 0.75. 
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